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Abstract. In this paper we give a detailed description of the random wavelet series representa- 
tion of real- valued linear fractional stable sheet introduced in By using this representation, 
in the case where the sample paths are continuous, an anisotropic uniform and quasi-optimal 
modulus of continuity of these paths is obtained as well as an upper bound for their behavior at 
infinity and around the coordinate axes. The Hausdorff dimensions of the range and graph of 
these stable random fields are then derived. 

1. Introduction and main results 

Let < a < 2 and H = (Hi, . . . , if at) S (0, 1)^ be given. We define an a-stable random field 
X = {X {t),t G R N } with values in R by 

Xo(t)= [ h H (t,s)Z a (ds), (1.1) 



where Z a is a strictly a-stable random measure on R with Lebesgue measure as its control 
measure and f3(s) as its skewness intensity. That is, for every Lebesgue measurable set A C R-^ 
with Lebesgue measure \n(A) < oo, Z a (A) is a strictly a-stable random variable with scale 
parameter Xj\r(A) 1 ^ a and skewness parameter (1/X^(A) ) f , 0(s)ds. If j3(s) = 0, then Z a is a 
symmetric a-stable random measure on R^. We refer to [20 . Chapter 3] for more information on 
stable random measures and their integrals. Also in (|1,1|) . 

N 1 1 

K (t, S) = K J] { (tt ~ Si)l l ~ ~ (-81)+'—} , (1.2) 
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where k > is a normalizing constant such that the scale parameter of Xq(1), denoted by 
ll-^o(l)IU) equals 1, t + = max{t,0} and 0° = 1. Observe that, if H\ = ■ ■ ■ = = i, Xq is the 
ordinary stable sheet studied in [9|] . In general, the random field Xq is called a linear fractional te- 
stable sheet defined on R^ (or (N, 1)-LFSS for brevity) in R with index H. LFSS is an extension 
of both linear fractional stable motion (LFSM), which corresponds to the case where N = 1, and 
ordinary fractional Brownian sheet (FBS) which corresponds to a = 2, that is, to replacing the 
stable measure in (11 . If) by a Gaussian random measure. 

We will also consider (N, <i)-LFSS, with d > 1, that is a linear fractional a-stable sheet defined 
on R N and taking its values in R rf . The (N, (f)-LFSS that we consider is the stable field X = 
{X(i),i G R N } defined by 

X(t) = (X 1 (t),...,X d (t)), VteR N , (1.3) 

where X\ , . . . , X& are d independent copies of Xq . It is easy to verify by using the representation 
(II. ip that X satisfies the following scaling property: For any N x N diagonal matrix A = (dij) 
with an = ai > for all 1 < i < N and = if i 7^ j, we have 

{X(At), t G R N } i {(ll«f ) X(t), * G K N ], (1-4) 

where = denotes the equality in the sense of finite dimensional distributions, provided that the 
skewness intensity satisfies (3 (As) = (3(s) for almost every s G R^. Relation (jl.4p means that the 
(N, ciY-LFSS X is an operator-self-similar [or operator-scaling] random field in the time variable 
(see (g, [23]). When the indices H\, . . . , Hn are not the same, X has different scaling behavior 



along different directions. This anisotropic nature of X makes it apotential model for various 
spatial objects, as is already the case for anisotropic Gaussian fields (Qj and 0]). We also mention 
that one can construct (N, d)-stable random fields which are self-similar in the space variables in 
the sense of QQ. This will not be discussed in this paper. 

Similarly to LFSM and FBS, see for instance fll. 13. fill. 11^. 21 . 24], there are close connections 
between sample path properties of LFSS and its parameters H and a. In this article we study 
some of these connections. In all the remainder of this paper we assume that the sample paths 
of X are continuous, i.e. min(i^i, . . . , i?j\r) > V a - For convenience we even assume that 

1/a < Hx < ■■■ < H N . (1.5) 

Of course, there is no loss of generality in the arbitrary ordering of H\, . . . , H^. Observe that, 
since H G (0, 1)^, condition (jl.5p implies a > 1. 

Let us now state our main results. Theorem Q] below is an improved version of Theorems 1.2 
and 1.3 in 0]. Relation (|1.6p provides a sharp upper bound for the uniform modulus of continuity 
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of LFSS, while Relation (jl.7p gives an upper bound for its asymptotic behavior at infinity and 
around the coordinate axes. 

Theorem 1. Let £Iq ^ e the event of probability 1 that will be introduced in Corollary^ Then for 
every compact set K, C R , all w G and any arbitrarily small rj > 0, one has 

sup — l*.(..uQ-*.(«.")l < „ (1 g) 

..^Ef.il^-'^'-^ti + liogl^-^ll) 2 ^' 

and 

SUP jr; ; : ; < CO. (1.7) 

tem nf=i \tj\ Hj (1 + I log md 1 /^ 
The following result can be viewed as an inverse of (|1.6p in Theorem [TJ 

Theorem 2. Ze£ 6e i/ie event of probability 1 t/iai ura'M be introduced in Lemma ["EH Then for 
all uj G fig, a// vectors u n 6 M^" 1 u>i£/i non-vanishing coordinates, any n = 1, . . . , N and any real 
numbers y\ < y2 and e > 0, one has 

\Xo(s n ,u n ,u)) - X (t 
sup _ , = 00, (1.8) 

s n ,t n e[vi,V2] |s n -t n | Jf "- 1 /«(l + |log|s n -t n ||) /a e 
where, for every real x n , we have set (x n ,u n ) = (u±, . . . , u n -i, ^n+i ; • • • , un)- 



Observe that Theorems Q] and [2] have already been obtained by Takashima 21] in the particular 
case of LFSM (i.e., N = 1). However, the proofs given by this author can hardly be adapted to 
LFSS. To establish the above theorems we introduce a wavelet series representation of Xq and 
use wavelet methods which are, more or less, inspired from [3|. It is also worth noticing that the 
event in Theorem [2] does not depend on u n . This is why the latter theorem cannot be obtained 
by simply using the fact that LFSS is an LFSM of Hurst parameter H n along the direction of the 
n-th axis. 

The next theorem gives the Hausdorff dimensions of the range 

X([0,l] N ) = {X(t): te[o,i] N } 

and the graph 

GrX{[0,l] N )={(t,X(t)):tG[0,l] N } 

of an (JV, d)-LFSS X. We refer to [10(] for the definition and basic properties of Hausdorff dimen- 
sion. 

The following result extends Theorem 4 in [2I to the linear fractional stable sheets. Unlike 
the fractional Brownian sheet case, we remark that the Hausdorff dimensions of -X"([0, 1]^) and 
GrX([0, 1]^) are not determined by the uniform Holder exponent of X on [0, 1]^. 
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Theorem 3. Let the assumption M.5\) hold. Then, with probability 1, 

dim H X([0,l] Ar )=min U ^— (1.9) 

I i=1 i ) 

and 

r k h N i } 

dim H GrX([0, 1]") = min j £ If + N ~ k + ( X " H ^ d > 1 ^ k ^ N ^Jf\ 

l ~ x l ~ l (1.10) 

w/tere 53°=i 35 : = 0- 

Remark 4. The second equality in (jl.lOj) can be verified by using (11. 5ft and some elementary 
computation; see 0]. 

In light of Theorem [3] it is a natural question to consider the Hausdorff dimensions of the 
image X(E) and graph GrX(E), where E is an arbitrary Borel set in WL N . As shown by Wu 



and Xiao 



22f | for fractional Brownian sheets, due to the anisotropic nature of X, the Hausdorff 



dimension of E and the index H alone are not enough to determine dim H X(.E). By combining 



the methods in Wu and Xiao 



22] and Xiao 



23f | with the moment argument in this paper we 
determine dim H X(£') for every nonrandom Borel set E C (0, oo) ; see Theorem 1211 

We end the Introduction with some notation. Throughout this paper, the underlying parameter 
spaces are R^, R+ = [0, oo)^ or "L N . A typical parameter, t G R^ is written as t = (t\, . . . , tjy) ° r 



t = (tj) whichever is more convenient. For any s, t £ R^ such that Sj < tj (j = 1, . . . , N), the set 



[ s > *] = T\f=i i s ji i s ca lled a closed interval (or a rectangle). Open or half-open intervals can be 
defined analogously. We will use capital letters C, C\, C2, . . . to denote positive and finite random 
variables and use c,c±,C2,... to denote unspecified positive and finite constants. Moreover, C 
and c may not be the same in each occurrence. 

Acknowledgment This paper was finished while the third author (Y. Xiao) was visiting the 
Statistical & Applied Mathematical Sciences Institute (SAMSI). He thanks Professor Jim Berger 
and the staff of SAMSI for their support and the good working conditions. 

2. Wavelet expansion of LFSS 

The goal of this section is to give a detailed description of the wavelet representations of LFSS 
Xq. First we need to introduce some notation that will be extensively used in all the sequel, 
(i) The real- valued function tp denotes a well chosen compactly supported Daubechies wavelet 
(see 0, [li3]). Contrary to the Gaussian case the fact that tp is compactly supported will 
play a crucial role in the proof of Theorem [2] (see the proof of Part (b) of Proposition [L 
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(ii) For any I = 1, . . . , N, the real- valued functions tfi Hi and ip~ Hl respectively denote the left- 
sided fractional primitive of order Hg + 1 — 1/a and the right-sided fractional derivative 
of order Hi + 1 — 1/a of tft, which are respectively defined for all x G R by 

= f(x- y)f- 1/a ^(y) dy and ^{x) = * f (y - x)% a - Hl ^{y) dy. (2.1) 

Observe that the functions ip He and tp~ He are well-defined, continuously differentiable 
and well-localized provided that tp has sufficiently many vanishing moments (and thus is 
smooth enough). By saying that a function : 1R — > ]R is well- localized we mean that 



sup(l + \x\) 2 {\cj){x)\ + \(j>'(x)\} < oo . (2.2) 



(iii) {e jk , (j,k) G Z N x Z N } will denote the sequence of random variables defined as 



N 



e ^= I f\h je/a ^ h Si-k e )\ Z a (ds). (2.3) 

They are strictly a-stable random variables all with the same scale parameter 

N/a 



hj,k\\a = { I \m\ a dt 

and skewness parameter 

Pj,k = Hk\\a a I TT {2 je i> <a> (V e Si ~ h)}P(s) ds 



N 



= 1 

where x <a> = |x| Q sgn(x) which is the number having the same sign as x and absolute 
value \x\ a . Moreover, if L > is a constant such that the support of tp is included in 
[-L, L], then for any integers p > 2L, any r G {0, . . . ,p — 1} N and j G Z JV , {ej jr +k p ; k G 
Z, N } is a sequence of independent random variables. 

A consequence of the above properties of the sequence {e^fc, (j, k) G Z N x Z N } is the following. 

Corollary 5. There exists an event of probability 1 such that, for any i] > 0, for all to G fig 

and all j,k € Z N x Z N , 

N 

|e ijfc M| <C(a;)n{(l + lid) 1/Q+,? (l + l^l) 1/Q log 1/a+, '(2 + |A ;i |)} , 
z=i 

where C is a finite positive random variable. 
Proof. We apply Lemma |2"51 □ 
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It is worth noticing that, for every I = 1, . . . , N, the functions ip He and ip~ He can be defined 
equivalently to (|2.ip . up to a multiplicative constant, but in the Fourier domain by (see e.g. fiol ] ) 

5*(0 = e'^^'-^^A (2.4) 

and 

= e^S^^-Va+l)! |^|^-l/a+l^_ ( 2 5 ) 

It follows from ParsevaFs Formula, (|2.4I) . (12. 5p and the orthonormality (in L 2 (R)) of the sequence 
{2-J'/ 2 V'(2 J ' • -k), j,k £ Z} that ip H * and satisfy, for all (J,K) £ Z 2 and (J', if) £ Z 2 , up to 

a multiplicative constant, 

/ *p Hl (2 J x - K)^- H z(2 J 'x - K>) dx = 2~ J 5(J, K; J 1 , K 1 ), (2.6) 
Jr 

where 6 (J, K; J' , K') = 1 when (J,K) = (J',K') and otherwise. By putting together (|2.4p . 
(|2.5p and the fact that = 0(£ 2 ) as |£| — > 0, another useful property is obtained: for every 
I = 1, . . . , N, the first moment of the functions ip Hl and ijj~ He vanish, namely one has 

ifj He (u)du = [ ifj- He (u)du = 0. (2.7) 

We are now in position to state the main results of this section. 

Proposition 6. Let fij be the event of probability 1 that will be introduced in Lemma UM. For 
every n £ N, M > and t £ R^ we set 

N 

U nM (t)= Yl ^ H) e jtk J{{^ Hl {2Hi-h)-^ Hl {-h)} , (2.8) 

where the random variables (j, fc) £ Z^ x Z N } are defined by A2.3\) and 

D% M = fc) £ Z N x Z w : for all 1 = 1,..., N < n and \k t \ < A/2 n+1 | . (2.9) 

Then for every uj £ £l\ the functional sequence (Z7 n ,Af( - 5 <^))neN is a Cauchy sequence in the Holder 
space C 7 (/C) /or every 7 £ [0, i?i — 1/a) and compact set K, C [— M, M]^. We denote its limit by 

N 

E a-^e^HW^tj - fo) - • 

(j,k)£Z N xZ N 1=1 

Proposition 7. With probability 1, the following holds for all t £ R^ 

N 

X (t) = E 2_ °' ,H>e ^n {V' H, (2 , '*i - **) - • ( 2 - 10 ) 
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Remark 8. By the definition of Xq and by Proposition [6l both sides of (|2,10p are continuous in 
t with probability 1. Hence, to prove Proposition it is sufficient to show that 

N 

E 2 ~ 0vH>e i,*II {^K&ti - h) - 4> H <(-k)} ,teR N 

K {j,k)eZ N xZ N 1=1 

is a modification of Xq. This is a natural extension of the wavelet series representations both 
of LFSM and FBS (see and will be called the random wavelet series representation of 

LFSS. 

Assume for a while that Proposition [6] holds and let us prove Proposition [7J 
Proof of Proposition [7J Let us fix I G {1, ... , N}. For any (ji, k{) G Z x Z and si G R we set 

^^(s^W^Si-fc,). (2.11) 

Since {ipj^k^j G Z, /c G Z} is an unconditional basis of L a (R) (see [3j) and, for every fixed t\ G R, 
the function ^ (t, - s/)+ i_1/a - (-s/)+ n ~ 1/Q G L a (R) H L 2 (R), one has 

(ti - ^)?" 1/Q - (-^)?~ 1/Q = EE KuttoWiiJnW . ( 2 - 12 ) 

where the convergence of the series in the RHS of (|2.12p , as a function of si , holds in L a (R) . Next 
by using the Holder inequality and the I? (R) orthonormality of the sequence 
jl G Z, ki G Z}, one can prove that 

Kl, jt , kl (tl) = 2^-^ [ {(n-st)^- 1 ^ -(-s^'-^^st-kOdsi 

JR 

= 2~ jl > H ' {i) Hl {2Hi-ki)-i) Hl {-h)} . (2.13) 

By inserting (|2.12p into (jl.ip for every I = 1, . . . , N, we get that for any fixed i G R^, the series 
([2.10p converges in probability to Xq (t) and Proposition [7J follows from Remark [8l □ 

From now on our goal will be to prove Proposition [6l We need some preliminary results. 
Proof of Proposition O For the sake of simplicity we suppose that N = 2. The proof for the 
general case is similar. The space C 7 (/C) is endowed with the norm 

||/|| 7 = sup|/(x)| + |/| 7 with |/| 7 = sup lf \* ] ~ f ^ )l , 
xe/C x^yeK. If - y\\ 7 

where || • || denotes the Euclidean norm in R 2 . For every n G N we set = (Z 2 x Z 2 ) \ L> 2 m . 
Let us define F n (x, y) = F n (x, y; ip Hl ; M, cp, 5, (3, rj) and E(x, y) = E(x, y; cfi, 5, (3, rf) by 



F n (x, y) = A n (x, y) + B n (x, y) , 
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where A n (x,y) and B n (x,y) are defined in Lemma [26] in the Appendix, and 



E(x,y)= y 2- ^ (2 X ~ K) -^ y - K)l (3 + \J\)V^(3 + \K\)V"+y. 

{J,K)& 2 1 y| 

Using ()2.8p . the triangle inequality and Lemma l22l one has for any n,p £ N and s\,s 2 ,ti i t 2 G 
[— M, M], denoting by Y\ the product over indices I = 1,2, 

|£4+p,A/(si, S2) — U n+Pi M(ti,t2) — U^m(si, S 2 ) + U ns M{ti,t 2 )\ 
(\si-t^ + \s 2 -t 2 \^ 

^ n - (jH) n ^ Hi (y Sl - h) - ^(-h)) - n (^'(^t, - *:,) - ^(-fc,)) 



(j'.*)G^ +p , M \ r, n,M 



< -\3,H) 



(\si-t 1 \ 2 + \s 2 -t 2 \ 2 )P/ 2 



(\si-h\ 2 + \s 2 -t 2 \^/ 2 



C C 



\si-h\P 

x \ip H2 (2 j2 s 2 - k 2 ) - ?p H2 (-k 2 )\ + similar term 
where C is a positive random variable. Observing that, for any non-negative array (oj,fc)(j,fc)ez 4 , 



£ a i> k - 
(i,fc)6A c l 



E E + E E 

(ji,fci)ez 2 YD„,M (i2,fc 2 )ez 2 (hM)& 2 O2, 



we thus get 

|^n+p,M(si, S2) — U n+Pi M(tl,t 2 ) — U n ^M(si,S 2 ) + J7 n ,Af *2)| 



< c 



(|si-tl| 2 + | S2 -t 2 | 2 )^ 2 

F n (si, ti; V^ 1 ; Hi, A r?)^(s 2 , 0; H 2 , 0, 77) + E{ Sl M^ Hl \ Hi, (3, r/)F n (s 2 , 0; ^ 2 ; # 2 , 0, r?) 



+ ^(s 2 ,t 2 ; i^ H2 ;H 2 ,f3, v)E(h,0; t/F 1 ; Hi, 0, rj) + £(s 2l t 2 ; *p H2 ;H 2 ,[3, rj)F n (ti,Q; ^ ; #1, 0, 77) 

By Lemma [26l we have that sup XiJ/g r MjM i F n (x,y) —> as n — ► 00 and sup x . i3/£ [_m,m] F(x,y) < 
00; hence the last display yields that sup„ >0 |{7 n +p,M — U U; M\-y — > as n — > 00. Observing that 
t^n,M vanishes on the axes, the same result holds with | • L replaced by || • || 7 and Proposition [6] 
is proved. □ 

Remark 9. Proposition [6] is much easier to prove in the Gaussian case. Indeed, in this case, 
using the fact that the e^'s are independent Af(0, 1) Gaussian random variables one can easily 
show that the sequence (Lfn,M}raeN is weakly relatively compact in the space C(/C). We refer to 
the proof of Proposition 3 in [3(] for more details. 
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From now on we will always identify the LFSS Xq with its random wavelet series representation 

3. Uniform modulus of continuity and behavior as \te\ — > or oo 

The goal of this section is to prove Theorem [TJ An immediate consequence of Proposition [6] 
is that Xq is locally C 7 for any 7 G (0, H\ — 1/a), almost surely. Theorem Q] completes this 
result by providing a sharper estimate on the uniform modulus of continuity, see (|1.6j) . and the 
behavior at infinity and around the axes, see (jl.7p . As in our note B, these s are obtained 
by using the representation (|2.10p . However, we improved the modulus of continuity estimate by 
relying on the independence present in the coefficients {e,- (j, k) G Z N x Z^}, see Lemma [23l If 
this independence is not taken into account, an alternative result (i.e., Lemma [2"2"j) may be used, 
resulting in a less precise estimate. The latter result holds in a quite general framework since 
they can be extended to a more general class of random wavelet series, see Remark 1101 below. 
Proof of Theorem Q3 It follows from (|2.10p . Corollary [5] and Lemma [271 that for every u £ S]q 
and every s,t G K, the triangle inequality implies 

N 

\X Q (s,iv) - X (t, u)\ < ^ |-^o (*i> • • • > tn-i, s n , . . . , sat; oj) — X (ti, . . . ,t n , s n+ i, . . . , sn;uj)\ 



n=l 

^H e ,i/ a , v {se;^ Hi 



jv / 7h — 1 \ / is 

< ciME n T ^W^;v^))x II T * 

n=l ^ e=i ' v= n +i 

xS HnA/a>ri (t n ,s n ;i) H ") (3.1) 

N 

< C 2 (U) £ K ~ S n r~ 1/a (1 + I log \t n - S n \\f /a+ "\ 

n=l 

This shows (fL6|) . 

Similarly, using (|2.10p , Corollary [5] and Lemma [271 we obtain, for every oj G VLq and every 
t G R, 

N N 

\X (t,u)\<C 3 (u)]jT He> y a „(tf,i; H *)<^ (3.2) 
The proof of Theorem Q] is finished. □ 

Remark 10. Clearly Proposition [6] holds more generally for any process Y = {Y(t), t G R^} 
having a wavelet series representation of the form 

N 

Y (t) = E c ^hk II {<t>i(V l ti - h) - H-h)} , 

(j,k)ez N xZ N 1=1 
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where the <^'s are well-localized functions, {cj k ,j,k G Z N } is a sequence of complex-valued 
coefficients satisfying \cj k\ < c2~^ >H ^ (c > being a constant) and {A^fe, j, k G 1< N } is a sequence 
of complex-valued random variables satisfying sup^ k E[|Aj ) / c | !/ ] < oo for all < v < a. We can 
also show that (jl.7p holds with probability 1 for such a process Y. In contrast, for this more 
general class of process, we cannot show (jl.6p but a less precise estimate for the uniform modulus 
of continuity. Namely, as announced in our note Q|, with probability 1, 

\X (s,u>) - X (t,u)\ 

SU P Ki 77 T, < 00 

S ,teKEf=ih--i;l^ 1/Q - r? 

for all compact sets K, C R-^. 

4. Optimality of the modulus of continuity estimate 

The goal of this section is to prove TheoremEJ For every n G {1, . . . , N} and (J n , k n ) G N x Z, 
let Gj ny k n = {Gj ny k„(u n ), u n G R^" 1 } be the a-stable field defined as the following wavelet 
transformation: 

G jn>kn {u n ) = 2^ l+H ^ [ X (s n ,u n )^~ H "(2^s n - k n ) ds n , (4.1) 

JR 

where the notation (s n , u n ) is introduced in Theorem [2j By using (|1.7p and the fact that the 
wavelet ip~ Hn is well- localized, the process {G? fc n (u), it G K^^ 1 } is well-defined and its tra- 
jectories are continuous, almost surely. The proof of Theorem [2] mainly relies on the following 
Lemmas [11] and [T2l 

Lemma 11. Let fig be the event of probability 1 in Corollary [5| and let n G {1, . . . , iV}. Suppose 
that there exist (u n ,u n ) G R^, p > 0, e > and uj G fig smc/i i/ioi 

|X (s n ,2 n ,a;) - X (t n ,2n, w)| 

SUp 1 -l/ a -t < °°' { > 

s n , *n£ [«„ -p,u n +p] \s n — t n |- f/ ' l_1 / a (l + | log | S n — t n 1 1) ^ 

T/ien one has 

limsup (j n 2- J '' l ) 1 / Q max{|G jn , fcn (S n ,t i ;)| : k n G Z, K - 2~^k n \ < p/8} = 0. (4.3) 

Lemma 12. Lei fij!j be the event of probability 1 defined as fij!j = fig n fi?i, where fig and fi?i are 
respectively the events defined in Corollary^ and Lemma[T3[ For all u G fi^, n G {1, . . . , iV}, a// 
integers j n G N, rea/ numbers z\ < z<i and all < t\ < T2, one has 

liminf (j n 2~^ a inf max \\G jn>kn (u n , w)|; fe n G [2***i, 2>'»z 2 ] n zl > 0. (4.4) 

jn->00 «nG[Tl,T 2 ] JV - 1 ^ J 

Before proving these lemmas, we show how they yield Theorem [2j 
Proof of Theorem [21 For the sake of simplicity we only consider the case where u n have 
positive and non-vanishing coordinates. The general case is similar. Suppose ad absurdum that 
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there exists u G such that (jl.8p is not satisfied. Then, for some n G {1, . . . , iV}, there exists 
u n G M^ -1 with positive and non vanishing coordinates, some real number u n , p > and e > 
arbitrary small such that (|4.2|) holds. By Lemma [TH this implies (|4.3|) . Then the conclusion of 
Lemma [121 leads to a contradiction. This proves Theorem [2l □ 
Proof of Lemma IX XL Let j n G N and k n G Z be such that 

\u n -2-i"k n \< p/8. (4.5) 

It follows from (|4.ip and (|2.7|) that Gj ni k n (u n ,ui) can be written as 

2i»(i+H») A ^ (« ni u n) t 1 )) - X (2-^*v,,«n, u))rl>- Hn C2? n s n - k n ) ds n . 
Hence, we have 

IG^^.W)! < 2^ 1+H ») / |X (a )l ,U rl ,w)-Xo(2-^fc n ,«n J w)| ^"^(^Sn-Ml^n 

= 2>^ 1+ ^){A in , fen (S n ,a;) + S in , fcn (2 n ,u;)} , (4.6) 

where 

A j n ,k n (u n ,oj) = / |X (s n ,2 n ,w) -X (2- J "fe„,S n ,cj)| |^- ff "(2 J "s n - A; n )|ds n (4.7) 

J\s„-u n \<p/2 

and 

Bj n ,k„(u n ,uj) = / |X (s n ,2 n ,w) -X (2- J '' l /c n ,S n ,u;)| |^-^"(2 J "s n -A: n )|ds n . (4.8) 

J |s„— u„|>p/2 

Let us now give a suitable upper bound for Aj n; k n (u n , w). It follows from (j4.7j) and (|4.2p that 
^j n ,fen(^,w) is at most 

CsM / |s n - 2-^fe„| H "~ 1/a (1 + I log|s n - 2-^k n \\y 1/a ~ e U~ Hn (2^s n - k n )\ds n . (4.9) 
Jr 

We claim that 

sup f \x\ Hn ~ 1/a {l/j n + I log 2 - (log|x|)/i n |)- 1/Q - e |^(x)| dx < 00 (4.10) 

and differ its proof after we have shown (|4.3p . 

By setting x = 2 Jn s n — k n in the integral in (|4.9p and using ()4.10p . one obtains, for all j n > 1 
and k n G Z satisfying (|4.5p . 

^(^w) < Q(a;)2^(- 1 -^ +1 /«) J --V«- . (4 . n) 
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In order to derive an upper bound for Bj n ^ n (u n , uj) , we use the fact that ip~ Hn is a well- localized 
function and f|4.5j) to get 

B j„,kn(un,u) < c I \X (s n ,u n , u) - X (2~ jn k n ,u n ,uj)\ (l + \2 jn s n - k n \) ds n 

J\s„-u„\>p/2 V 7 



< c / \X (s n ,u n ,u) - Xq(2 jn k n ,u n ,uj)\ 

J \s n —Un\>p/2 



S„—U n \>p/2 

x(l + 2 jn (\s n - u n \ - \u n - 2~ jn k n \^j 2 ds n 
< c2~ 2jn \X (s n ,u n ,uj) - X (2~ jn k n ,u n ,u;)\\s n - u n \~ 2 ds n . 

J\s n -u n \>p/2 

This last inequality, together with (jl.7p . implies that, since ui G fig, 

B jn>kn (u n ,u;) < C 7 (cu)2~ 2 ^, 

where Cj is a random variable that does not depend on the integers j n and k n satisfying (|4.5p . 
Hence, putting together the last inequality, (|4.1ip and f|4.6[) one obtains (|4.3|) , 

Finally, to conclude the proof of the lemma, it remains to show (|4.10p . We separate the 
integral in (|4.1U|) into two domains, \x\ > 2 Jn / 2 and \x\ < 2 Jn / 2 . We bound (l/j n + | log 2 — 
(log \x\)/j n |) 1 6 from above by jn a+IL on the first domain, and by ((log2)/2) e on the 
second domain, yielding that the integral in (|4.10p is at most 

Jn /a+e [ \x\ H "- 1/a \^- H -(x)\dx+ ((log2)/2)~ 1/cM f \x\ Hn ~ 1/a \$- H »(x)\dx . 

Using that H n — 1/a G (0, 1) and that tp~ Hn is well localized, we thus get (|4.10p . □ 

In order to prove Lemma [T2l we first prove a weaker result, namely the following lemma. 

Lemma 13. There exists fi?ij an event of probability 1, such that for all uj G fi?i, n G {1, . . . ,N} 
and real numbers M > 1, z\ < z 2 , < t\ < r 2 , one has 

liminf (j^-^^uinjn, M; Zl , z 2 ;t 1 ,t 2 ;uj) > , (4.12) 

where 

v(n,j n ;M;zi,z 2 ;Ti,T 2 ;uj) 

(4.13) 



min max <^ \G jn , kn (M Jn k n ,uj)\; k n G 

k n e[M^ Tlj M^ T2 ] N - 1 nz N - 1 



2 ]n Zl , 2^z 2 



n; 



In order to prove Lemma [Till we need to show that the random variables Gj n ^ n (u n ) satisfy 
some nice properties, namely the following proposition. 

Proposition 14. Let u n G R^" 1 be an arbitrary fixed vector with non-vanishing coordinates. 
Then the following results hold: 
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(a) {Gj nj k n (u n ), (jnj^n) £ N x Z} is a sequence of strictly a-stable random variables with 
identical scale parameters given by 



\Gjn,k n { u n)\\a 



L a (R) IJ II 



\H t -l/a 
> + 



xHi-l/a 



(4.14) 



(b) Let L > be a constant such that the support of tp is included in [—L,L]. Then for all 
integers p > 2L and j n > 0, {Gj niqnP (u n ); q n S Z} is a sequence of independent random 
variables. 



Proposition O is in fact a straightforward consequence of the following proposition and the 
fact that any two functions s n t— > tp(2 3n s n — q n p) with different values of q n have disjoint supports. 

Proposition 15. For every vector u n with non-vanishing coordinates and for every (j n ,k n ) 6 
NxZ one has almost surely 



G jn ,k n @n) = [ \2^(2^s n - k n ) II - si)+~ 1/a - dZ a (s) . (4 



.15) 



Proof of Proposition 1151 As in (|3.2p . we have 



sup 



nf=i \tj\ H *{i+ 1 log i^i 



< oo. 



(4.16) 



It follows from Propositions [7] and EJ (|4.16p . the Dominated Convergence Theorem, (I2.6p . (12.7 
and (12.13P that for any u n S IR^ -1 one has almost surely, for any increasing sequence (D. 



'm/m£N 



of finite sets in Z x Z such that U m D r 



G Jn ,k n (un)= lim 2 J " 



-Un,H n ) 



N 



f nU H '(2 j 's/-^)-^(-^)]^(2 J "Sn-A: n )d Sn (4.17) 



= m lhn oo^ n«^.*«( U ') £ ( 3 nJn);(*n,*n)' 

where «| 7, is defined in (|2.13p . On the other hand, it follows from (|2.12p that 



Si 



)?~ 1/Q }= e n^^nvw^, (4.i8) 



1 = 1 
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where for all fixed u n G M. N 1 the convergence of the series in the RHS (|4.18f) . as a function of 
s G R N , holds in L a (R N ). Next using (|4TT8l) and §Z3§ one has, for every fixed u n G M^" 1 , 



dZ a {s) 



(i„,fcn)ez 2 ( JV - 1 ) 



where the convergence of the series holds in probability. Finally, putting together (|4.17p . (|4.19p 
and (|2.1ip . one obtains the proposition. □ 

We are now in position to prove Lemma [T3l 
Proof of Lemma fT3l For any constants M, c\ > 0, n G {1, . . . , N}, integ 

jn ^ and rational 

numbers r\ < r 2 , < 0\ < 02 and C > 0, let F(n,j n ) = T(n,j n ; M, c\\ r 2 ; #i, 02] C) be the event 
defined as 



T(n, j n ;M, a; ri,r 2 ; 0i,02) = {oj : v{n, j n ; M; n, r 2 ; 0i, 2 ; u) 
First we will show that, there exists c\ large enough such that 

w(r(n,j n ;M, Cl ;r 1 ,r 2 ;6 1 ,9 2 )^ < 



oo. 



l,...,N t>l (j„,fc„)eNxZ UnGlOi,^]^- 1 



for ui 7^ 0, we have 

C2 := min inf inf 

n- 

Observe finally that 

v(n,j n ;M;ri,r 2 ;0i,02;uj) 

>^ min max\\G jniqnP (M~ jn k n ,uj)\; q n e 

k n £[Mine 1 ,M Jn d 2 ] N - 1 nZ N - 1 { 

It follows from Proposition Q3] and (|4.22p and this inequality that 

r(n, Jn ))<E n 



2 -j„^-l/a| 


(4.20) 




(4.21) 


; with \ui\ and 


non-zero 


> i) > . 


(4.22) 


-2 Jn n 2- 7n r 2 " 


nzj. 












(4.23) 



< C3 M( JV - 1 ^(l-c 2 j n 2-^/ci) 



where the summation is taken over all k n G [M-?"^, JVP"^]^ 1 H 1 and the constants c 2 , C3 
and C4 do not depend on . Using the last inequality one can prove that (|4 . 2 1 1) holds by choosing 
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c\ > large enough. Hence the Borel-Cantelli Lemma implies that, for such a constant c±, 



'( U Pi r c (n,j n ;M,c 1 ;r 1 ,r 2 ;9 1 ,0 2 )\ =1, 



where r c (n, j n ; M, c±, r±, r 2 ; 6*1, # 2 ) denotes the complement event of T(n, j n ; M, c\\ r±, r 2 ; 0\, 62)- 
But this implies that the event 

lo : ]hrLmf(j n 2-i n ) 1 / a v(n,j n -,M;r 1 ,r 2 -,6 1 ,e 2 ;u) > 

j n ->oo 

has probability 1. Finally setting S7 2 as the intersection of such sets over |(M; r%, r 2 ; 0%, 62) G 

Q 5 : M > 0, r\ < r 2 and < 9\ < # 2 |, one obtains the lemma. □ 

The following proposition will allow us to derive Lemma [12] starting from Lemma [131 Roughly 
speaking it means that the increments of the random field {Gj n ^ n (un), u n G [ti,^]^ -1 } can be 
bound uniformly in the indices j n and k n . 

Proposition 16. Let £Iq be the event of probability 1 that was introduced in Corollary 0. Then 
for any reals z\ < z 2 , < n < r 2 and r] > arbitrarily small, there exists an almost surely 
finite random variable C$ > such that for every n G {1, . . . , N}, j n G N, k n G [2 Jn zi, 2 Jri 2 2 ] J 
S n G [ri, t 2 ] w n G [ri, t^^" 1 and a; G Qq, one has 

\G jn , kn (u n ,u) - G jn , kn (v n ,u)\ < C 8 ( w ) 2^ K ~ vi\ H '- 1/a - v - (4.24) 

Proof. Lemma [23 applied to (|3.1|) shows that, for all oj G £1^ and any r\ > 0, there exists 
C(u) > such that, for every n G {1, ... , iV}, s n G R, u n G [n, T2] n ~ 1 and u n G [ti, r 2 ] iV_1 , 



|X ( Sn ,2 n , W )-X ( Sn ,^,u;)| <C(cu)( ^l^-^l^-^-M l^l^^l + llogd^DI) 1 ^. (4.25) 

Let £ > be arbitrary small and consider the integral 

I{j n ,k n ) = 2>» I (1 + \s n \) H ^ \^- H -{2^s n - k n )\ ds n . 

By setting x = 2 Jn s n — k n we derive that 

sup max I(jn,k n ) = sup max / I 1 + 2 Jn \x + k n \ I \ip n (x) \dx 

j n e®k„.e[Z>nzi,2jnz 2 ] j n m k n e[2inzt,2in Z2 ] J R \ J 

f / \ H n +C . „ 

< / 1 + |x| + max{|z!|, |z 2 |} \if)~ Hn (x) \dx < 00. (4.26) 

The inequality p~24j) then follows from ([44]), ([4351) and ([4~26l) . □ 

We are now in position to prove Lemma [T2l 
Proof of Lemma 1121 We set 

inf max \\G jn , kn (u n ,u)\; k n G [2 3 "2i, 2^z 2 ] n zl. (4.27) 
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In view of Lemma [13] it is sufficient to show that there exists 7 > small enough and M > 
such that, for all n E {1, . . . , N}, u 6 and reals z\ < z 2 , < t\ < T2, one has 



lim 2- J " (1/a - 7) 



u(n, j n ;M;z 1 ,z 2 ;T 1 ,T 2 ;uj) - v(n, j n ] zi, z 2 ;t 1 ,t 2 ;u) =0. (4.28) 



As the function fj n (-) = max | |Gj nj fc n (-, w)|; & n £ [2 Jn ;zi, 2 Jn Z2]| is continuous, there exists 
v%,(jn) G [n, t^' 1 such that 

/i„(^0n)) = inf {f jn (u n ); u n G [ti.ts]^ 1 }. (4.29) 

Moreover, when j n is big enough, one has for some k®(j n ) € [M Jn ri, M J "T2] Ar ~ 1 H Z^ -1 , 

||M">4°(j n ) - ^(in)||oo < M-f» . (4.30) 

Then it follows from Proposition \W\ that there exists a constant C5 > (independent of (j n ,k n )) 
such that the following inequality holds 

\G jntkn (M-^k° n (j n ),u;) - G jnjkn (vP n (j n ),u)\ < c 5 

The last inequality implies that 

/,„(M-^^(j n )) < / in (u°(j„)) + c 5 ^F^-'M. (4.31) 

By using (|4.29|) and (|4.31j) one obtains that 

f jn (u° n (jn)) < mm {f jn (M-^k n ); k n G [M**n, M^]"" 1 } 

Let us choose M large enough so that 

H N -l/a log M 
fli - 1/a < log 2 ' 

and then, using (jl.5|) . we choose 77 > and 7 > small enough so that 2^ H "M-^ Hl ~ 1 / a ~^ = 
o(2~ Jn ( 1//Q! ~ 7 )) as j n — > 00. Finally combining this with (|4.3ip . we obtain f|4.28[> . This proves 
Lemma [T2l □ 

5. Proof of Theorem [3] 

As usual, the proof of Theorem[3]is divided into proving the upper and lower bounds separately. 
The proofs of the lower bounds rely on the standard capacity argument and the following Lemma 
[TT1 However, the proofs of the upper bounds are significantly different from that in [3J, due to 
the fact that both dim H X([0, 1]^) and dim H GrX([0, 1]^) are not determined by the exponent 
'or the uniform modulus of continuity of X. Our argument is based on the moment method in 



(4.32) 



12|. Combining this argument with the methods in [23J, we are able to determine the Hausdorff 
dimension of the image X(E) for all nonrandom Borel sets E C (0, 00) . 
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We start by proving some results on the scale parameters of the increments of real-valued LFSS 
Xq between two points (i.e., Xq(s) — X^(t)) and over intervals; see Lemmas 1171 and 1181 below. 
Combining the latter with the maximal moment inequality due to Moricz [l^j, we derive sharp 
upper bounds for the moments of the supremum of Xq. 

Lemma [T71 is an extension of Lemma 3.4 in Q| for fractional Brownian sheets. Since d(s,t) := 
||Xo(s) — Ao(t)|| a can be used as a pseudometric for characterizing the regularity properties of 
Xq via metric entropy methods (cf. 2JJ, Chapter 12]), these results will be useful for studying 
other properties of LFSS X as well. 



Lemma 17. For any constant e > 0, there exist positive and finite constants cq and C7 such that 
for all s,t£ [e, 1] N , 



N 



N 



C6^\s£-t £ \ ' <\\Xq(s) - X (t)\\ a < C 7 ^2\s £ -t e \ l . 

1=1 e=i 



(5.1) 



Proof. To prove the upper bound in (15. lj) . we use induction on N. When N = 1, Xq is an 
(H, a)-linear fractional stable motion and one can verify directly that (|5.ip holds as an equality. 
Suppose the upper bound in (15. ip holds for any linear fractional stable sheet with n parameters. 
We now show that it holds for a linear fractional stable sheet Xq with n + 1 parameters. 

It follows from the representation (jl.ip that, for any s,t £ [e, l] n+1 , ||Xo(s) — Xo(t)||° is a 
constant multiple of the following integral: 



„ n+l 1 1 n+1 

jUn {<*-•*)?-' -<-dr s } -n 

1=1 £=1 



/ \ Hp 

{st-n) + l a 



-re) 



dr 



< c 



n _ n 

n fa - - c-^)? - -} -ufa- ^) 



(-re) 



dr 



|(Wl - r n+l) + 



(-r n+ i)^" +1 "| 



n+l 



+ c 



dr 



< c 



1 1 
H{(^-r,)?~-(-r^~} 

t=i 

x / \(t n+1 - r n+ i)+" +1 a - (s n+1 - r n+ i] 
JR 1 

n \ a ) 

L,\st-tt\ Ht ) + |Wl - S n +l\ Hn+ia \ , 



+1 



where, in deriving the last inequality, we have used the induction hypothesis, the fact that the 
function t 1— > f R {{t — r)^_ — {—r)^_ 1//Q } a <ir is locally uniformly bounded for any H > 1/a 
and that, by a change of variable r n+ \ = t n+ \ + \t n+ \ — s n+ i\u, the last integral in the previous 
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display is less than \t n+ x — s n+ i\ aHn+1 up to a multiplicative constant. Hence we have proved the 
upper bound in (|5.ip . 

For proving the lower bound in (|5.1|) . we define the stable field Y = {Y(t),t G M^} by 



Y(t)= / h H (t,r)Z a (dr), (5.2) 
J[o,t] 

where the function h H (t,r) is defined in (11. 2j) . Then by using (II. ip again we can write 

\\X ( S )-X (t)\\>\\Y(t)-Y(s)\\. (5.3) 



To proceed, we use the same argument as in [3|, pp. 428-429] to decompose Y as a sum of 
independent stable random fields. For every t £ [e,l] N , we decompose the rectangle [0, t] into the 
following disjoint union: 

N 

[0, t] = [0,e] N U \jR(tj)UA(e,t), ( 5 - 4 ) 

i=i 

where R(tj) ={r£ [0, 1]^ : < rj < e if z 7^ j, e < 7j < i,} and A(e, t) can be written as a union 
of 2^ - N- 1 sub-rectangles of [0,<]. It follows from (JO} and fl53D that for every t £ [e, 1]^, 

TV 

Y(t) = [ h H (t,r) Z a (dr) +V / h H {t,r) Z a {dr) + ( h H (t,r) Z a (dr) 

J[0,e]» ~{JR(t 3 ) JA(e,t) 

N 

:= Y(e,t) + ^2Y j (t) + Z(e,t). (5.5) 

Since the processes {Y(e,t), t £ R N }, {Yj(t), t £ R N } (1 < j < N) and {Z(e,t), t £ R N } are 
defined by the stochastic integrals with respect to Z a over disjoint sets, they are independent. 
Only the lj(i)'s will be useful for proving the lower bound in (|5.1|) . 

Now let s,t G [s, 1] N and j £ {1, ...,N} be fixed. Without loss of generality, we assume 
Sj < tj. Then 

||^(*)-^WHa= / {h H (t,r)-h H (s,r)) a dr+ [ h a H (t,r)dr, (5.6) 
JR(sj) JR(s h tj) 

where R(sj,tj) = {r G [0, 1]^ : < n < e if i / j, Sj < Tj < tj}. By {ESD and some elementary 
calculations we derive 

\\Yj(t)-Yj(s)f a > / h a H (t,r)dr 
JR{s :j ,t j ) 

= [ Uitk-rk)^"- 1 f\t 3 -r 3 r H ^dr (5-7) 



Mi 

~-^ > C I tj Sj I ^ ■ 
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where c > is a constant depending on e, a and (1 < k < N) only. The lower bound in (|5.ip 
follows from (|51fll and ([521) • □ 

Our next lemma determines the scalar parameter of the increment of Xq over any interval 
[ s 't\ = IT^=i [ s ji~tj]- Recall that the increment of Xq over [s,t], denoted by Xo([s,t\), is defined 
as 

Xo([s,t}):= Yl (-if^Xodsj + Sjitj-sj))). (5.8) 

5e{0,lp 

This corresponds to the measure of the set [s, t] by interpreting Xq as a signed measure defined 
by X ([0,t\) = X (t) for all t G R N (convention: [0,tj] := [tj,0] if tj < 0). It may be helpful to 
note that for N = 2, we have 

X ([s,t}) = X (t) - X (( Sl ,t 2 )) - X ((t 1 ,s 2 )) + X (s). 

Similarly, we will denote the increment of the function hff(-,r) over [s,t] by hu([s, t], r). 

Lemma 18. For any interval [s,t] = YljLi i s j^j]' we have 

N 

\\Xo([s,t])\\a = U^-^ aHj - ( 5 - 9 ) 
Proof. Since the kernel hff(t, s) in (jl.2p is a tensor product, it can be verified that 
ll*o(MC= I \h H ([s,t],r)\ a dr 



N i i 



air 



(5.10) 



v I=i 1 

N 

11(0 *,•)«"'••. 

This proves Lemma [THJ □ 

In order to estimate E[sup teT \Xn(i)— Xg(a)\] for all intervals T = [a, b] C [e, 1]^, we will make 
use of a general moment inequality of Moricz [17[ for the maximum partial sums of multi-indexed 
random variables. This approach has the advantage that it is applicable to non-stable random 
fields as well. Another way for proving Lemma [20l below is to establish sharp upper bounds for 
the tail probability P[sup teT \Xo(t) — Xp (a) | > u] by modifying the arguments in [18|. 

First we adapt some notation from [17|] to our setting. Let k € N N } be a sequence of 
random variables. For any m £ Z+ (Z+ is the set of nonnegative integers) and k £ N^, let 
R = R(m, k) = (to, to + k] D Z+ , which will also be called a rectangle in Z+, and we denote 

5(i?) = 5(m, fc) = and M(i?) = max |S(m,g)|. (5.11) 



20 ANTOINE AYACHE, FRANQOIS ROUEFF, AND YIMIN XIAO 

It can be verified that M(R) < maxgcR p(Q)| < 2 N M(R), where the maximum is taken over all 
rectangles Q C R. Let f{R) be a nonnegative function of the rectangle R with left-lower vertex 
in . We call / superadditive if for every rectangle R = R(m, k) the inequality 

/(%) + f(R j2 ) < f(R) (5.12) 

holds for every 1 < j < N and 1 < qj < fcj, where 

Rji = i?((mi,.. . ,mjv), (fci, ■ ■ • , Qj, . . .,fejv)) 

and 

-Rj2 = R((mi, . . .,mj-i,mj + qj,m j+ i, . . .,m N ), (fci, . . . fcj - . . . , fejv))- 

In other words, Rji U i?j2 = R is a disjoint decomposition of i? by a hyperplane which is perpen- 
dicular to the jth axis. Together with the nonnegativity of /, (|5.12p implies that, for every fixed 
m G Z5: , f(R(m, k)) is nondecreasing in each variable fcj (1 < j < iV). 



The following moment inequality for the maximum M(R) follows from Corollary 1 in 17]. 



Lemma 19. Let (3 > 1 and 7 > 1 be given constants. If there exists a nonnegative and su- 
peradditive function f(R) of the rectangle R in such that E[|5'(i?)| 7 ] < f^(R) for every R, 
then 

v.[M{Rp] < (l) N {i - 2^~^) Nl AR) ( 5 - 13 ) 



for every rectangle R in 



7N 



It is useful to notice that the constant in f)5 . 13|) is independent of R. Applying Lemma [19] to 
the linear fractional stable sheets, we obtain 

Lemma 20. Let the assumption tl.5\) hold. Then there exists a positive and finite constant c$, 
independent of the skewness intensity /3(s), such that for all rectangles T = [a, b] C [e, 1]^, 

e( sup \X (t) - X (a)\) < c 8 V(6j - aj) H '. (5.14) 
Proof. We prove this lemma by using induction on N. In the case of N = 1 it is well known 



(cf. [12j or 2l|) that (|5.14p holds. Observe that the term 5 = (1) in the sum appearing in (|5.8p 
is Xo(t), and since J2se{o i} N (~^) N ~^ tSi = 0' we have 

Xo([s,t})= (-if-^iMisj + Sjih-sjm-Xois)) 

5&{0,1} N 

= X (t) - X Q (s) + Yl (- 1 ) 7V " El5l {^o((s J +^(^-s i )))-X ( S )} . (5.15) 
5e{o,i} Ar \{i> 
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For every 5 G {0, 1}^ \ (1), there is some n G {1, . . . , N} such that 5 n = 0- Observing that, using 
the notation introduced in Theorem[2] u n \— > Ao(a n , u n ) is an [N — 1, 1)-LFSS as defined by (jl.ip 
and (|1.2|) but with N replaced by N — 1, k multiplied by a constant only depending on a n and 
bounded independently of a n (since a n G [e, 1]) and a modified skewness intensity (3. Hence using 
the induction hypothesis, we have, for every 5 G {0, 1}^ \ (1), 

Ef S up|X ((a i + ^(t i -a i )))-X (a)| N ) < c £> y - a,-)*' . (5.16) 

Applying (15.150 with s = a and (|5.16j) . the bound (I5.14|) is implied by 

N 

E(sup |X ([a,t])|) < c Hibj - aj ) H i , (5.17) 

teT 3=1 

which we are now going to prove. This is where Lemmas HH1 and [19] will be applied. 
For all n G N we define a grid in [a, b] with mesh 2~ n by the collection of points 

r n (p) = {aj+ Pj (bj - aj )2- n ), p G R(0, (2 n )) = {1, . . . ,2*}" . 

For each p G -R(0, (2 n )), we define the random variable ^ p to be the increment of Xq over the 
elementary rectangle with upper-right vertex T n (p), [r n {p — (l)),T n (p)]. Interpreting X as a 
signed measure, we get that for any rectangle R(m, k) C R(0, (2 n )) for 

^ ^ = X ([r n (m), r n (m + *)]). (5.18) 

We are now ready to prove (|5.17|) . By the continuity of the sample function Xo(i) and the 
monotone convergence theorem, since the set U n >i{r n (p) : p G R(0, (2 n ))} is dense in [a, b], it is 
sufficient to show that for all integers n > 1, 

e( max |X ([a, r n (p)])| > ) < c 9 TT(^ " "if** ( 5 - 19 ) 

where eg > is a finite constant independent of [a, 6] C [e, 1]^ and ra. 

It follows from Lemma [23] in the Appendix that for any strictly a-stable random variable Z 
with scale parameter 1 and every < 7 < a, we have E(|Z| 7 ) < cio, where c\q depends on a 
and 7 only. This fact, (|5.18|) and Lemma [T5] imply that for every 1 < 7 < a and every rectangle 
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R = R(m,k) C i?(0,(2 n )), 



E 



E 



X ([T n (m), T n (m + A;)]) 



< cio 



< 



JV 

n 

N 



kj(bj aj / 



Hi 



Hj/m-iHn 



(5.20) 



where c n = c{^\ 



For every rectangle -R = R(m, k) included in i?(0, (2 n )), let 



AT 



kjibj-aj) 



2" 



Hj/Hi 



Note that, under assumption (jl.5j> . we have a S (1, 2), i?ia > 1 and Hj > Hi for j = 1, . . . , N. 
Hence the inequality x Hj ^ Hl +y H ^ Hl < (x+y) H ^ Hl for all x, y > implies that / is superadditive. 
We take 7 G (1, a) such that ft = 7 -Hi > 1 and apply Lemma [T9l to derive 



E 



sup 

fcei?(o,{2™)) 



pG-R(0,fe) 



< C12 



A? 



nfe-«i) ,/ ' /,/i 



TV 



Hi 7 



(5.21) 



II /' 

where C12 > is a finite constant independent of [a, b] and n. This proves (|5.17p and thus Lemma 
M □ 

We now proceed to prove Theorem [3J 
Proof of Theorem [3J We only prove (|1.9|) , which is done by modifying the proof of Theorem 4 
in Q| and by making use of Lemmas [T7] and [20j The formula (jl.lOp can be proven using similar 
arguments and we leave it to the interested reader. 

First we prove the lower bound in (|1.9[) . Let e G (0, 1) be given and let I = [e, 1] N . We will 
prove that for every < 7 < mm{d, Yle=i TT e }' dim H X(7) > 7 almost surely. By Frostman's 
theorem (see e.g. [101 ] pages 64 and 65), it is sufficient to show that we have 

1 



E 



dsdt < 00, 



(5.22) 



Ijjj \\X(s)-X( 

where || • || denotes the Euclidean norm in M. d . 

It is known that for any ci-dimensional distribution function F in M, d with characteristic function 
<p and any 7 > 0, we have 



27 /2-l r /7 

V2 



p p+00 p 

/ \\x\\-^F(dx) = (2vr)~ d / 2 / u^du 
jR d Jo Jm 



exp 



ip(ux) dx. (5.23) 
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This equality can be verified by replacing 99 in the right side of (|5,23p by its expression as a 
Fourier integral and then performing a routine calculation. Applying (|5.23j) to the distribution 
of the stable random variable £ = (X(s) — X(t)) /\\X(s) — X(t)\\ a and using Fubini's theorem, 
we obtain 

E(||eH" 7 ) < cm / exp ( - dx r u 7 " 1 exp ( - ci B M a ||a;H du 

= Ci6 / exp I ||x|| 1 dx < 00, 

jR d V 2 / 

where the last integral is convergent because 7 < d. Combining (|5.24p with Lemma [T71 yields 

1 

Ijjj \\X(s) -X(t)p " Jjjj (£f =1 \ at - U\^y 



E / / 7rF7 -^ v ,^^, dsdt < j j ^ N -7777^7 dsdt < °°> ( 5 - 25 ) 



where the finiteness of the last integral is proved in 0, p. 432]. This proves (|5.22p and hence the 
lower bound in (jl.9p . 

To prove the upper bound in (jl.9p . we use the covering argument in [3] and Q]. Since clearly 
dim H X([0, 1]^) < d a.s. and Hausdorff dimension is cr-stable irj], it is sufficient to show that for 
every e S (0, 1), 



N 

j=l 1 



N 1 

d\m H X([e,l] N ) <EtT a ' S - (5 - 26) 



This will be done by using a covering argument. 

For any integer n > 2, we divide [e, 1]^ into m n sub-rectangles {R n ,i} with sides parallel to 
the axes and side- lengths n~ l l Hi (j = 1, . . . , N), respectively. Then 

m n < cn 3=1 h j (5.27) 

and X(Je, 1] N ) can be covered by X{R n ^) (1 < i < m n ). Denote the lower-left vertex of R n ^ by 
a n ^. Note that the image X{R n ^) is contained in a rectangle in M. d with sides parallel to the axes 
and side lengths at most 2sup sgRn . |X&(s) — X k (a n ^)\ (k = 1, ... ,d), respectively. Hence each 
X{Rn,i) can be covered by at most 

"2sup sei?n . \X k (s) - X k (a n! i)\ 



n 

k=l 



n 1 



+ 1 



cubes of side-lengths n 1 . In this way, we have obtained a (\^dn ^-covering for X([e, 1] N )- 
By Lemma [201 we derive that for every 1 < i < m n and 1 < k < d, 

E( sup \X k (s) - X k (a n ,i)\) <cn~ l . (5.28) 



21 



ANTOINE AYACHE, FRANQOIS ROUEFF, AND YIMIN XIAO 



It follows from (15.27L (I5.28|) and the independence of X±, . . . ,Xd that for any 7 > Ylj=i 7F' 



we have 



« En 

i=i k=i 



2 sup sei?n . \X k (s) - X k (a n!i 



-1 



+ 1 



n 



(5.29) 



< cn in' 



as n 



00. 



This and Fatou's lemma imply that dim H X([e, 1]^) < 7 almost surely. By letting 7 j X^fci TT 
along rational numbers, we derive (|5.26p . This completes the proof of Theorem [3j □ 

The above method can be extended to determine the Hausdorff dimension of the ima ge X (E) 
for every nonrandom Borel set E C (0, 00)^, thus extending the results in Wu and Xiao 22( and 



Xiao 



23j for anisotropic Gaussian random fields to (N, (i)-LFSS. 



For this purpose, let us first recall from [23J the definition of a Hausdorff-type dimension which 
is more convenient to capture the anisotropic nature of X. 

For a fixed (Hi, . . . , H jy) £ (0, 1)^, let p be the metric on defined by 

N 

p(s,t) =^2\sj -tj\ H] , Vs,teR N . (5.30) 



For any > and E C R , define the /^-dimensional Hausdorff measure [in the metric p] of E 
by 



f OO OO N 

H$(E) = liminf I ^(2r n )^ ■■ Ec[jB p ( r n ), r n < 8 >, 

^ n=l n=l J 



(5.31) 



where B p (r) denotes a closed (or open) ball of radius r in the metric space (R , p). Then Hp 
is a metric outer measure and all Borel sets are Hp -measurable. The corresponding Hausdorff 



dimension of E is defined by 



dim££ = inf {(3 > : H P p (E) = 0}. 



(5.32) 



We refer to 



23l ] for more information on the history and basic properties of Hp and dim£ . 
Theorem 21. Let the assumption i f 1.5]) hold. Then, for every nonrandom Borel set E C (0, oo) N , 

dim H X(£) = min {d; dim^E} a.s. (5.33) 



Proof. The proof is a modification of those of Theorem [3] above and Theorem 6.11 in [231 ] . For 
any 7 > dim^E, there is a covering {B p (r n ),n > 1} of E such that Y^=i(^ r nV ^ 1- Note that 
X(E) C U^ =1 X(Bp(r n )) and we can cover each X(B p (r n )^j as in the proof of Theorem [3j The 
same argument shows that dim H X(E) < 7 almost surely, which yields the desired upper bound 
for dim H X(E). 
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By using the Frostman lemma for Tip (Lemma 6.10 in 23J) and the capacity argument in the 
proof of Theorem [3l one can show dim H X(E) > minjd; dim^E 1 } almost surely. We omit the 
details. □ 

Appendix A. Technical lemmas 

The following lemma allows to control the growth of an arbitrary sequence of strictly a-stable 
random variables having the same scale parameter. 

Lemma 22. Let {e\, A £ Z, d } be an arbitrary sequence of strictly a-stable random variables 
having the same scale parameter. Then, there exists an event of probability 1, such that for 
any r] > and any us £ fi*, 

d 

|€ A H|<CHn(3 + |A,|) 1/0 ^, (A.l) 
i=i 

where C > is an almost surely finite random variable, only depending on r\. 

Proof This lemma simply follows from the fact that for any v £ ((1/a + f?) -1 , ol) one has 

E ( sup —, ~~~~ ~7~ ) < c Y TT(3 + \\A)- u W a+ ri < oo. 

□ 

Lemma 23. Let a £ (0,2). There exists a constant cyj depending only on a such that for 
any strictly a-stable random variable Z with scale parameter \\Z\\ a > and skewness parameter 
(3 £ [—1, 1] and all t > ||^||q, 

c^ 1 \\Z\\% t~ a < P(|Z| > t) < ci7 \\Z\\ a a t~ a . (A.2) 

Let N > 1. Suppose now that {Zj^, ji > 1, kg > 2 for £ = 1, . . . , N} is a sequence of strictly 
a-stable random variables such that 

(i) For all j £ (N \ {0})^, {Zj )k , h > 2 for i = 1, . . . , N} are independent; 

(ii) For all j £ (N \ {0})^ and k £ (N \ {0, 1})*, \\Z jjk \\ a < 1. 
Then, with probability 1, one has, for any 7 > 0, 

sup Ylj; 1/a ' 7 k; 1/a log- 1 /^ ki : u > 1, h > 2 for I = 1,... < 00 . (A.3) 



Proof Relation (|A.2|) follows from Property 1.2.15 in [20J. Let us now show (|A.3|) for N = 1, the 



proof for iV > 1 is similar. By using (|A.2[) . we obtain, for all j > 1 and n > 1, 
P(max{|%,|, . . . , \Z j>n \} > Uj , n ) < 1 - (1 - c l7 uj«) n , 
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where Uj iTl = j l / a+ i n l / a log 1//Q+7 n. Defining n m = [exp(m)], we obtain 



E 



EE 1 . 



ma,x{\Zj >2 \,---,\Zj }nm \}>Uj tn7 
j>l m>l 



E P ( maX {l%2|, • • • , \Zj,n m \} > U j)Um j < 00 • 
j>l m>l 



Thus the random variable £\>i E m >i 1 max{|z J , 2 |,...,|z J , nm |}> % , nm is a.s. finite. As a consequence 
there exists an a.s. finite positive random variable C such that 

max{|% 2 |, • • • , \Zj,nJ\} < Cu jinm for all j > 1, m > 1 . 

Let m(fc) be the unique integer satisfying n m nA < k < n m ^ +l . Thus for all j > 1, k > 2, we 
have 

fel < ^> m(fe)+1 = ^i Va+7 log 1/a+7 (n OT(fc)+1 ) , 

Observe now that we have, for all k > 2, 

n m ( fc ) + i < exp(m(/c) + 1) < e (n m(fc) + 1) < e (k + 1) . 

Relation (|A.3[) follows from the last two displays. □ 

Lemma 24. For any 7 E [0, 1) and n > 0, i/iere exists a constant c > smc/i t/iai, /or a// 
^(2 + |u - k\)~ 2 (l + |fc|) 7 log"(2 + |fe|) < c(l + |u|) 7 log"(2 + |«|) . 

fcGZ 

Proof Put k' = [u] — k, where [u] is the integer part of u. Hence 

E (1+ f2 i r| 1 ° g -il^ |fc|) = E( 2 + I- - M + k'\)-Hl + \[u] - k'W log^(2 + \[u] - k'\) 

< ^(l + |fc , |)- 2 (2 + |n| + |A : , |) 7 ^(2 + ^1 + 1^1) . 

The result then follows by observing that (2+|u| + |/c'|) 7 < (l + |n|) 7 (2 + |fc'|) 7 , log"(2+M + |A/|) < 
c log" (2 + |u|) log" (2 + \k'\)} and 7 - 2 < -1. □ 

Lemma 25. Let 9 ^ and 7 6 K. Sef c := En>o 2 "' e|n ( 1 + n )' 7 ' < 00 ■ T/ien / or 

any no < ni 

in {0,±1,±2, ... ,±00}, 

A fl f2 nofl (l + I n 1) 7 if6<0 

y 2 ne (l + |n|) 7 <c{ V ' W (A.4) 
nt^ (^(l + lml) 7 i/fl>0. 

Proof. Take e.g. 8 < and write 

£ 2^(1 + |n|) 7 < 2^(1 + |n |r E ^ (^T^rY " 

n=n m>0 ^ ' °' ' 

Now observe that 

1 l + \m + n \ 

< ■ — ■ — < 1 + m 



1 + \m\ 1 + |n | 
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so that sup no E m >o 1 rnd ( i TT|Sr) 7 < °° for any 7 G R ' D 

Lemma 26. For any M > 0, n > small enough, 6 G + n, 1), /? £ [0, <5 — 1/a — n), any 
well-localized function (ft and x, y G R, /ei A n (s, y) := A n (s, n; M, 0, <5, /3, 77) be the quantity defined 
as 

A n (x,y)= Y Y 2~ JS ~ K) ~ ^ y - K)l (3 + I J|) 1 /»-H>( 3 + |if j) 1 /^ ( A .5) 

\J\<n\K\>M2 n + 1 1 1 



and Zei B n (x,y) := B n (x, y; (ft; 5, (3, rj) be the quantity defined as 

w 



4-0 (l , v) _ E l^'-fj^-^ p+W)^. (A. ( 



B„(x,»)= Y. E 2- J< l0(2Jl ~ tfi " K>l (3 + I J|)'^ + -(3 + |if |)'/«+-, (A.6) 

tuii/i toe convention that A n (x,x) = B n (x,x) = for any x £ R. These quantities converge to 0, 
uniformly in x,y G [— M, M], as n goes to infinity. 

Proof. Let x,y £ [—M,M] and Jo > — log 2 (2M) be the unique integer such that 

2'- 10 - 1 <\x-y\<2- Jo . (A.7) 

Let us first prove that A n (x,y) converges to 0, uniformly in x,y as n goes to infinity. From now 
on we suppose that J is an arbitrary integer satisfying \J\ < n. We need to derive suitable upper 
bounds for the quantity 

)J^. _ Tf \ — ■> 

\x — y\P 

\K\>M2 n + 1 1 y ' 

For this purpose, we consider two cases J < Jo and J > Jq + 1 separately. First we suppose that 

J < J . (A.9) 
Using the Mean Value Theorem, ([2T2]) . (|A.7|) and (|A.9j) one obtains that 

\<ft(2 J x -K)- <p(2 J y - K)\ < c2 J \x - y\ sup(3 + \u\)~ 2 

uei 

< c2 J \x-y\(2 + \2 J x- K\)~ 2 , 

where / denotes the compact interval with end-points 2 J x — K and 2 J y — K, whose length is at 
most 1 by (|A.7p and (|A.9|) . Next the last inequality and (|A.8P entail that 

A^,y)<c^-yf-e E iTW^W- <A ' I0 » 

\K\>M2 n + 1 y 1 U 

On the other hand, using that |x| < M and \ J\ < n, for all \K\ > M2 n+ , one gets 

<3 + |*|)V«« < (3 + jgD'/'+i ^ < e ( 1 + |y|)-(»-'/<-i). (A . U) 



(2 + |2 J x - K|) 2 ~ (2 + - M2 n ) 2 
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Putting together (|XT0]) . (TA"TLj) and (TATF]), one obtains that 

A { n J \x,y) < c2 J o(/3-i) 2 ^-«(i-i/"-r ? )_ (A. 12) 

Let us now study the second case where 

J > J + 1. (A.13) 
It follows from (lA?7j) . (fA~T3l) and (fO|) that 

4 J )(x,y) < 2 J ^ ^ {|0(2 J x - + \cf>(2 J x - K)\}(3 + \K\) 1 ^. (A.14) 

|A'|>M2™+ 1 

On the other hand, using (|2.2p and the fact that | J| < n one has, for any real u £ [—M,M] and 
any iTsZ satisfying > M2 n+1 , 

\cp{2 J u- K)\ <c{3 + \2 J u-K\)- 2 < c(3 + \K\ - M2 n )~ 2 < c 18 (3 + \K\)~ 2 . (A.15) 

Combining (|A.14p with (|A. 15[) one gets that 

A^\x,y) < ci9 2 J ^- n ( 1 - 1 / a -") . (A.16) 

It follows from COl) . (fOj) . (lATHl) and (TAlBl) that 



< c2 



-ra(l— l/a— rj) 



Jo oo 

2 Jo(/3-i) g 2 J ( 1 -^(3 + |J|) 1/a+r '+ 2 J ( /3 - <5 )(3 + |J|) 1 / Q+r ' 

J= — OO J=J0 + 1 

< c 2- n ( 1 - 1 /«-^) 2 J o(/ 3 - 5 )(3 + \J \) l l a+r i 

< C 20 2 -"( 1 - 1 /"-'?) ; 

where we used Lemma [231 to bound the series and then the fact 2~ Jo < 2M (see ()A.7j) ) . Since 
C20 does not depend on (x,y), the last inequality proves that A n (x, y) converges to 0, uniformly 
in x, y S [-M, M] as n goes to infinity. 

Let us now prove that B n {x,y) converges to 0, uniformly in x,y as n goes to infinity. In all 
the sequel J denotes an arbitrary integer satisfying \J\ > n + 1. First, we derive a suitable upper 
bound for the quantity 

)J ry. _ it\ _ sf,(2 J ° 

w 



Bi ) _ 1^-/0-0(2^-/01 + /a+ „ 



X — 7/ 



As above, we distinguish two cases: J < Jq and J > Jo + 1. First we suppose that (|A.9|) is 
verified. As in (jAlOj) . we have B { n J \x,y) < c2 J \x - y] 1 ' 13 Ea g z( 3 + |# |) 1/a+r? (2 + |2 J x - K|) -2 . 
Next, using (|A.7|) and Lemma [241 and the fact that \x\ < M, one obtains that 

^ J )(x,y) < c2 J+Jo ( /3 - 1 )(l + 2 j ) 1/q+,? - (A.18) 
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Now let us suppose that (|A.13p is verified. By using this relation, (|A.7j) . the triangle inequality, 
(|2.2p . Lemma [2^1 and the fact that x, y £ [—M,M], one gets 



4 J )(x,y) < 2 J ^{\ < P(2 J x-K)\ + \ ( / ) (2 J y-K)\}(3 + \K\) 1 / a + r ' 

< c2 J/3 {( 3 + l 2 ^ - K \)~ 2 + ( 3 + l 2 ^ - -^l)~ 2 }( 3 + \K\) 1/a+v 



< c2 J/3 {(l + 2 J \x\) 1/a+r > + (1 + 2 J |y|) 1 / a+r '} 

< C 2 j( .f 3+1 / a +v) . (A.19) 
Since 2~ J ° < M, for all n > log 2 (2M), we have — n < Jo, and thus, by (|A.18j) . 

^ 2- J<5 (3 + |J|) 1 / Q +'?B(- / )( a;)2/ ) < C 2 J °( /3 - 1 ) 2 J{1_5) (3+ |J|) 1/q+?? 

< c2 n ( 5 - 1 )(l + n) 1 / Q+ '' , (A.20) 
where we used Lemma [25] and 2 _J ° < M. Applying Lemma [25] with (|A.18j) and (|A.19j) yields 

^2" J5 (3+ \ J\) 1/a+ri Bi J \x,y) < c2 Jo ^ +1 / a+r i-V(3+\J \) 1/a+r > , (A.21) 

and for any n > Jo, 

^2" J,5 (3 + \J\) 1/a+r 'B^ J \x,y) < c 2 n ^ +1 / a+r i- 5 \3 + n) 1 / a+ri . (A.22) 

J>n 

Since /3 + 1/a + n - 5 < 0, the function t i-> 2 t ^ +1 / a+r '-^ (3 + t) 1 /"-^ i s decreasing for t large 
enough, and hence for n large enough, either n > Jq and we may apply (lA~22l) . or n < Jo 
and we may apply (|A.2ip whose right-hand side is smaller than the right-hand side of (|A.22j) . 
Hence (|A.22|) holds for all n large enough independently of Jo- This, with (|A.20p . shows that 
B n (x, y) converges uniformly in x, y, as n goes to infinity. □ 

Lemma 27. Let (f> be a well-localized function i.e. a function satisfying the condition For 
any 5 £ (0, 1), 7 £ (0, 5) and n > 0, define 

S S , 7 , v (x, y;<j))= £ 2- J5 \^(2 J x - K) - <p(2 J y -K)\(3 + | J|)^(3 + \K\)i log^(2 + \K\) 

(A.23) 

and 

T J)7)IJ (x;0)= £ 2" J5 |^(2 J x-K)-0(-K)|(3 + |J|)^'(3 + |K|)T'log^(2 + |^|). (A.24) 

(J,ir)ez 2 
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Then, there exists a constant c > 0, only depending on 5, 7 and eft, such that the inequalities 

Ss,~f, v (x, y\<t>) <c\y - x\ 5 ~^ [\y - x| 7 + |x| 7 + |?/| 7 ] 

x (1 + I log \y - x\ { \og 1+r, {2 + \x\)+ log 7+r '(2 + \y\)} (A.25) 

and 

T S)1)V (x;4>) < c(l + jloglxlD^lx) 5 (A.26) 
hold for all x, y £ R (with the convention that a x log 6 = for all a, b > 0). 

Proof. We only prove (|A,25p . the proof of (|A.26P is similar. By (|2.2p . there is a constant c > 
such that, for all J,K 6 Z and x, y G R, 

|</>(2 J x -if) - c/>(2 J ?/ - K)| < c {(2 + |2 J x - K|)" 2 + (2 + \2 J y - K\)~ 2 } . (A.27) 

The quantity \cp(2 J x - K) - cj)(2 J y - K)\ can be bounded more sharply when the condition 
2 J \x — y\ < 1 holds, namely by using (|2.2p and the Mean Value Theorem one obtains that 

\(p(2 J x - K) - (p(2 J y - K)\ < c2 J \x - y\ sup(3 + \2 J u - K\)~ 2 

< c2 J \x-y\(2 + \2 J x- K\)~ 2 , (A.28) 

where I denotes the compact interval whose end-points are x and y. From now on we will assume 
that x 7^ y (Relation (|A.25|) is trivial otherwise) and let Jo € Z be the unique integer satisfying 

1/2 < 2 Jo \y-x\ < 1. (A. 29) 

The inequalities ()A.27P and (|A.28|) entail that 

Ss,y, v {x,y,(/)) < c(A Jo \x-y\ +B Jo ), (A.30) 

where 

A Jo = ^22 J ^- s \2 + \2 J x-K\y 2 (3 + \J\y +r, (3 + \K\yiog^(2 + \K\) 

J<J K&L 

and 

B ->o = E E 2- J5 {(2 + |2 J x-if|)- 2 + (2 + |2 J y-Kr 2 }(3 + |J|)^(3 + |^|) 7 log^(2 + |^|). 

J> J KGZ 

Lemma l24l and Lemma [25] yield 

A Jo < c2 J °M (1 + |x| 7 2 J ° 7 ) (1 + |J |) 27+2?? log 7+ "(2 + |x|) 
and, since 7 — 5 < 0, 

Bj < c2- Jo8 (1 + (|x| 7 + \yp)2 J ^) (1 + I Jol) 2 ^ 2 " {log 7+ "(2 + |x|) + log 7+? ?(2 + \y\)} . 
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Inserting these two bounds into (|A.30j) and using (|A.29p . we get ()A.25|) and the proof is finished. 
□ 
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